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Abstract 

We consider the parameter estimation problem for the non-ergodic fractional Ornstein- 
Uhlenbeck process defined as dX t = 9X t dt + dB t , t > 0, with a parameter 9 > 0, where B 
is a fractional Brownian motion of Hurst index H £ (|, 1). We study the consistency and 
the asymptotic distributions of the least squares estimator 6t of 9 based on the observation 
{X s , s € [0,t]} as t oo. 
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1 Introduction 

We consider the Ornstein-Uhlenbeck process X = {X t , t > 0} given by the following linear stochas- 
tic differential equation 

X Q = 0; dX t = 9X t dt + dB t , t > 0, (1) 

where B is a fractional Brownian motion of Hurst index H > \ and 6 g (—oo, oo) is an unknown 
parameter. An interesting problem is to estimate the parameter 9 when one observes the whole 
trajectory of X. First, let us recall some results in the case when B is a standard Brownian motion. 
In this special case, the parameter estimation for 6 has been well studied by using the classical 
maximum likelihood method or by using the trajectory fitting method. If 9 < (ergodic case), the 
maximum likelihood estimator (MLE) of 8 is asymptotically normal (see Liptser and Shiryaev [9] , 
Kutoyants [8J). If 6 > (non-ergodic case), the MLE of 9 is asymptotically Cauchy (see Basawa 
and Scott [3], Dietz and Kutoyants [4J). Recently, in a more general context, several authors 
extended this study to some generalizations of Ornstein-Uhlenbeck process driven by Brownian 
motion (for instance, Barczy and Pap |2J). Similar properties of the asymptotic behaviour of MLE 
has also been obtained with respect to the trajectory fitting estimators (see Dietz and Kutoyants 

a)- 



When B is replaced by an a-stable Levy motion in the equation ([T]) , Hu and Long [B] discussed 
the parameter estimation of 9 in both the ergodic and the non-ergodic cases. They used the 
trajectory fitting method combined with the weighted least squares technique. 

Now, let us consider a parameter estimation problem of the parameter 9 for the fractional 
Ornstein-Uhlenbeck process X of (JXJ) . 

In the case 9 < (corresponding to the ergodic case), Hu and Nualart [7] studied the parameter 
estimation for 9 by using the least squares estimator (LSE) defined as 

a jlx s dx 

folds' ~ 

This LSE is obtained by the least squares technique, that is, 9 t (formally) minimizes 



+ ex. 



2 



ds. 



To obtain the consistency of the LSE 9t, the authors of [7] are forced to consider J Q X s dX s as 
a Skorohod integral rather than Young integral in the definition ([2]). Assuming J*X s dX s is a 
Skorohod integral and 9 < 0, they proved the strong consistence of 9t if H > s, and that the LSE 
9 t of 9 is asymptotically normal if H £ [A, §). Their proof of the central limit theorem is based 
on the fourth moment theorem of Nualart and Peccati [12] . 

In this paper, our purpose is to study the non-ergodic case corresponding to 9 > 0. More 
precisely, we shall estimate 9 by the LSE 9t defined in ([2]), where in our case, the integral J Q X s dX s 

is interpreted as a Young integral. Indeed in that case, we have 9t = „ r t * 2 -, which converges 

2 Jo X s ds 

almost surely to 9, as t tends to infinity (see Theorem [TJ. Moreover, it turned out that the path- 
wise approach is the preferred way to simulate numerically an estimator 9t ■ Our technics used in 
this work are inspired from the recent paper by Es-Sebaiy and Nourdin [5]. 

The organization of our paper is as follows. Section 2 contains the presentation of the basic 
tools that we will need throughout the paper: fractional Brownian motion, Malliavin derivative, 
Skorohod integral, Young integral and the link between Young and Skorohod integrals. The aim 
of Section 3 is twofold. Firstly, we prove when H > | the strong consistence of the LSE 9 tl that 
is, 9t converges almost surely to 9, as t goes to infinity. Secondly, we investigate the asymptotic 
distribution of our estimator 9t in the case H > \. We obtain that (see Theorem [5]) 

e et (d t -9 s ) ^ 29C{1) as t — > oo, 

with C(l) the standard Cauchy distribution with the probability density function ^jz^r; > x E R. 



2 Preliminaries 

In this section we describe some basic facts on the stochastic calculus with respect to a fractional 
Brownian motion. For more complete presentation on the subject, see [11], [T] and |10| . 
The fractional Brownian motion (B t ,t > 0) with Hurst parameter H £ (0,1), is defined as a 
centered Gaussian process starting from zero with covariance 

R„(t, s) = E(B t B.) = \ (t 2H + s 2H |t - sf H ) . 

We assume that B is defined on a complete probability space (CI, J 7 , P) such that T is the sigma- 
field generated by B. By Kolmogorov's continuity criterion and the fact 

E{B t -B s f = \s~t\ 2H - s,t> 0, 
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we deduce that B has Holder continuous paths of order H — e, for all e £ (0, iJ). 

Fix a time interval [0,T]. We denote by T-L the canonical Hilbert space associated to the 
fractional Brownian motion B. That is, Ti is the closure of the linear span £ generated by the 
indicator functions l[o,t] 5 t £ [0, T] with respect to the scalar product 

(l[o,t] 5 l[o,s]> = s). 

The application ip £ £ — > B(ip) is an isometry from £ to the Gaussian space generated by B and 
it can be extended to Ti. 

If H > i the elements of "H may be not functions but distributions of negative order (see |13j). 
Therefore, it is of interest to know significant subspaces of functions contained in it. 
Let \H\ be the set of measurable functions <p on [0,T] such that 



< oo. 



IMIftt, :=H(2H-l) / / ^(uJII^JHu-ul^-'dttdw 
Jo Jo 

Note that, if ip, ip £ \U\, 

E(B(cp)B(ip)) = H(2H - 1) / / <p(u)^(u)|u- w| 2H_2 dudw 



It follows actually from [T3] that the space \H\ is a Banach space for the norm and it is 

included in W. In fact, 

L 2 ([0,T]) C LTt{[Q,T\) C \H\ C H. 

Let Cg°(K n ,M) be the class of infinitely differentiable functions / : R" — > R such that / and 
all its partial derivatives are bounded. We denote by S the class of smooth cylindrical random 
variables F of the form 

F = f(B( n ),...,B(<p n )), (3) 
where n > 1, / e Cg°(R n ,R) and <pi,...,<p n £ H. 

The derivative operator D of a smooth and cylindrical random variable F of the form ([3]) is defined 
as the "H-valued random variable 

N Q, 

DtF = E q^.WW-), ...,B{!p n ))<pi(t) 

In this way the derivative DF is an element of L 2 (£l; T-L). We denote by D 1 ' 2 the closure of S with 
respect to the norm defined by 

ll^lli, 2 -^(ll^l| 2 )+^(ll^ll«)- 

The divergence operator 8 is the adjoint of the derivative operator D. Concretely, a random 
variable u £ L 2 (£l;H) belongs to the domain of the divergence operator DomS if 

E\{DF,u) u \<c\\F\\ L 2 (n) 

for every F £ S. In this case S(u) is given by the duality relationship 

E(FS(u)) =E(DF,u) n 

for any F £ D 1,2 . We will make use of the notation 

cT 



3( u ) = / u s5B S: u £ Dom(5). 
Jo 



In particular, for h £ H, B(h) = 6(h) = L h e SB, 
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For every n > 1, let be the nth Wiener chaos of B , that is, the closed linear subspace 
of L 2 (il) generated by the random variables {H n (B(h)), h £ H, \\h\\-n = 1} where H n is the nth 
Hcrmite polynomial. The mapping 7„(/i®") = n\H n (B(h)) provides a linear isometry between the 
symmetric tensor product H " (equipped with the modified norm || - H^O" = -^=||.||%gm) and H n . 
For every /, g £ H Gn the following multiplication formula holds 

E(I n (f)U9)) = n\(f )9 )n®«- 

Finally, It is well-known that L 2 (il) can be decomposed into the infinite orthogonal sum of the 
spaces H n - That is, any square integrable random variable F £ L 2 (fl) admits the following chaotic 
expansion 



F = E(F) + J2 Wn) 



n=l 

where the /„ £ are uniquely determined by F. 

Fix T > 0. Let /, g : [0,T] — > K are Holder continuous functions of orders a £ (0,1) and 
j8 £ (0, 1) with a + j3 > 1. Young P3] proved that the Riemann-Stieltjes integral (so-called Young 
integral) J Q f s dg s exists. Moreover, if a = ft £ and <j> '■ ^ — > i is a function of class 

C 1 , the integrals J Q gu)df u and J §^(/«, g u )dg u exist in the Young sense and the following 

change of variables formula holds: 



f ^{fu,9u)df u + [ ^(f u ,9u)dg u , < t < T. (4) 
Jo of J dg 



ft ax pt 

<P(fu9t) = <P(fo,9o) + 



As a consequence, if H > | and (u t , t g [0, T}) be a process with Holder paths of order a > 1 — H, 
the integral f Q T u s dB s is well-defined as Young integral. Suppose moreover that for any t £ [0,T], 
u t £ D 1 ' 2 , and 



^ J \D s u t \\t- s\ 2H ~ 2 dsdt <oc}j = 1. 



P 

Then, by [T], u € Bomb and for every i G [0, T], 

/ u s dB s = [ u s SB s + H(2H - 1) f [ D s u r \s - r\ 2H ~ 2 drds. (5) 
Jo Jo Jo Jo 

In particular, when ip is a non-random Holder continuous function of order a > 1 — H , we obtain 

/ p s dB s = [ cp s 6B 8 = B(<p). (6) 
Jo Jo 

In addition, for all <p, ip £ \H\, 



E 



[ ip s dB s [ ip s dB s J = H(2H - 1) / / ^"^Kcfo. (7) 

v jo Jo y Jo Jo 



3 Asymptotic behavior of the least squares estimator 

Throughout this paper we assume H £ (|, 1) and 8 > 0. Let us consider the equation ([T]) driven 
by a fractional Brownian motion B with Hurst parameter and 9 is the unknown parameter to 
be estimated from the observation X. The linear equation ([1]) has the following explicit solution: 

X t = e et f e- 9s dB sl t > 0, (8) 



o 
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where the integral J Q e dB s is a Young integral. 
Let us introduce the following process 



(9) 



^ ■= / e~ vs dB s , t>0. 
Jo 

By using the equation ([T]) and (j8]) we can write the LSE 9t defined in ([2]) as follows 

fix*da folds' 

3.1 Consistency of the estimator LSE 

The following theorem proves the strong consistency of the LSE 9 t . 

Theorem 1 Assume H g (^,1), then 

9 t — >9 almost surely 

as t — > oo. 

For the proof of Theorem Q] we need the following two lemmas. 

Lemma 2 Suppose that H > h. Then 

i) For all e G (0, H), the process £ admits a modification with (H — s)-Holder continuous paths, 
still denoted £ in the sequel. 

ii) — > S,oo '■= Jq°° e~ 0r dB r almost surely and in L 2 (VL) as t — > oo. 
Lemma 3 Let H > A. Then, as t — > oo, 

e~ 2et f X 2 s ds = e~ 2et f e 26s e s ds — > %■ almost surely. 
Jo Jo 20 

Proof of Lemma [2j We prove the point i). We have, for every < s < t, 

/ft v 2 

\ 2 _ T?l „-0r r 



t r t 

2H-2, 



Efo-cr = E [j e dB 

= H(2H-l)j j e- 8u e- ev \u~v 
< H(2H-1)J J \u-v\ 



dudv 

2H - 2 dudv 

E{B t -B s f = \t-sf" 



Thus, by applying the Kolmogorov-Centsov theorem to the centered gaussian process £ we deduce 
i). 

Concerning the second point ii), we first notice that the integral £oo = J e 8r dB r is well defined. 
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In fact, 



poo poo 

H(2H-1) / e' 9r e- es \r - s\ 2H - 2 drds 
Jo Jo 

poo 1>S 

2H(2H-1) dse~ 9s dre' dr \s - r) 2ff ~ 2 
Jo Jo 



u u 2H-2 



pOO pS 

= 2H{2H - 1) / dse- Ws / due 
Jo Jo 

POO />oo 

= 2H(2H-1) due 8u u 2H - 2 dse 

Jo Ju 
H(2H - 1) f°° _ 2H-2, 

" e u du 



-29s 



V r(2g - 1} = < °°' (10) 



with r denotes the classical Gamma function. Moreover, ^ converges to ^ in L 2 (f2). Indeed, 

/OO pC 
It 



poo />oo 

e e " " I r — s 



"' • ' K • ^ 2H - 2 drds 



2H-2 



/OO />S 
dse ~ 9S J t dre- 9r {s-r) 

poo ps — t 

= 2H{2H - 1) / d ue - 29s / due 6u u 2H - 2 
Jt Jo 

pOO pV 

= 2H(2H - l)e- 29t / dve- 26v / due 9u u 2H - 2 
Jo Jo 

pOO pOO 

= 2H(2H - l)e- 2et / due 6u u 2H - 2 / dve~ 2ev 

Jo Ju 
H(2H — T) 2gt f -0u„.2H-2, 



e leu 



l du 



mm _ 2et 

Q2H 

— > as t —> oo. 

Now, let us show that £ t — > ^ almost surely as t — > oo. By using Borel-Cantelli lemma, it is 
sufficient to prove that, for any e > 



sup f e- 0s dB s 

n>0 V"<t<n+1 Jt 



>e) < oo. (11) 



For this purpose, let i < a < H . As in the proof of [Theorem 4, pQ], we can write for every i > 

,-8s jn „— 1 / jo „— 6>s / / j_/„ _\ — a/„ j.\a:— 1 



2 

oo 



e- as dB s = C - L j dB s e-° s yj dr(s - r)~ Q (r - t) 

with c Q = J. (s — r)~ Q (r — t) a ~ 1 dr — /3(a, 1 — a), where /3 is the Beta function. 
By Fubini's stochastic theorem (see for example [11]). we have 

/oo />oo / />oo 

e~ Bs dB s = c- 1 J dr(r - t)^ 1 (j dB s e~ 9s (s - r)~ Q 



fi 



Cauchy-Schwarz's inequality implies that 
2 

e-°'dB K 



< c: 



jf° O (r _ t )2(a-l) e -0(r-t) dr ^ dre -e(r-t) jf 



C ^r(2a -l) e _ 20t r, re -^) 



02a-l 



/OO />OG 



dB,e- l "(« - r)- a e' 

2 



(s-r)- a e- 0(5 - r) 



Thus. 



sup 

n<t<n+l 



e~ &s dB s 



< 

On the other hand 



c- 2 r(2a-l)„- 29n „ fl , x 



Z 1 OO />00 

/ dre- ^-^ / dB s (s-r)- a e- 8 < s - r ) 



{^J\s-r)- a e- e ^ s - r ^dB s 
= H{2H-\) / - r )- a e' e{v '^ / du(u - r) 

Jr »/ r 

pOO pOO 

= H{2H-\)\ dvv- a e- Bv duu- a e- 6u \u-v\ 
Jo Jo 

/>00 f-V 

= 2H(2H-l) dvv- a e~ ev duu~ a e- Bu {v - u) 2H - 2 
Jo Jo 



a e- 9 ^\u-v\ 2H - 



2H-2 



2H(2H-1) dvv- a e-" v du{v-u)- a e 
Jo Jo 



a p -8(v-u) u 2H-2 



poo PV 

< 2H{2H-1) dvv- a e- 0v du{v-u)- a u 
Jo Jo 



a n ,2H-2 



duu 2tl -\l-u)- a 



= 2H(2H-1) y°° dtw 2ff - 2a - 1 e- fl " J 

= 2H(2H-l) g2g-2a := Ci{a,H,0) < oo. 

Combining this with the fact that e~ e ^ r ~ n ^dr = |, we obtain 



_E sup 

V n<t<n+l 



e~ &s dB< 



< C 2 (a,H,6)e 



-20n 



with 



C 2 {a,H,6) = 



c- 2 T(2a-l)e e 

Q 2a 



Consequently. 



n>0 



sup 

\n<t<n+l 



J e- 0s dB s >e^j 



> e) ' e 



n>0 



a,H, 6). 






/•OO 


sup 




\ n<t<n+l 





< e- 2 C 2 (a,H,9)J2 e ~ 29n < °°- 



n>0 



This finishes the proof of the claim (fTTI) , and thus the proof of Lemma [21 ■ 
Proof of Lemma [3j Using ([TU]) , we have 

Hence ^ ^ A/"(0, ^J^J 1 ^ ) and this implies that 

P(£oo = 0) = 0. (12) 
The continuity of £ entails that, for every t > 

* /•* f / N 



e 2Ws Cds > / e^'ffda > -e wt inf £ almost surely. (13) 
Furthermore, the continuity of £ and the point ii) of Lemma [5] yield 

Una inf £ I = almost surely. 
Combining this last convergence with (fT3j) and (fT2l) , we deduce that 

lim / e (£ds — oo almost surely. 
i-* 00 Jo 

Hence, we can use L'Hospital's rule and we obtain 

, e & 

lim — — ^ = lim — = — - almost surely. 

t-Hx> e 29t t^oo 29 29 

This completes the proof of Lemma [3] ■ 

Proof of Theorem [TJ Using the change of variable formula ((4]) , we conclude that 



2 

Hence 



1 e m g=9 f e 2es e s ds+ f e 9s ^dB s 



o jo 



J* e 2f)s e s ds 2e~ 2et £ e 29s £ 2 ds 
Combining this with Lemma [2] and Lemma [3j we deduce that 9t — > 9 almost surely as t — > oo. ■ 

3.2 Asymptotic distribution of the estimator LSE 

This paragraph is devoted to the investigation of asymptotic distribution of the LSE 9t of 9. We 
start with the following lemma. 

Lemma 4 Suppose that H > i. Then, for every t > ; we have 

dB s e 0s / dB r e- 6r = / dB s e 6s / dB r e- 0r - / SB s e- 8s / 6B r e 9r 
(i Jo Jo Jo Jo Jo 



H(2H - 1) f dse- 9s [ dre 
Jo Jo 



0rt s _,2H-2_ 
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Proof. Let t > 0. By the change of variable formula (HJ) 

dB,e 8s / dB r e~ 9r = / dB,e 6s / dB r e~ 6r - I dB, e - 6s I dB r e 0r . 



On the other hand, according to ([5]) and ((6j), 

dB s e- 9s I dB r e 9r 

i 



f SB s e- 9s f 5B r e 6r + H(2H ~l) f dse' es f dre e 
Jo Jo Jo Jo 



9r\„ „{2H-2 



\s-r\ 



which completes the proof. ■ 

Theorem 5 Let H > i be fixed. Then, as t — > oo, 

e et fg t _ g \ i^ 2 0C(l), 



with C(l) the standard Cauchy distribution. 

In order to prove Theorem [3] we need the following two lemmas. 

Lemma 6 Fix H > h. Let F be any a {B} -measurable random variable such that P(F < oo) = 1. 
Then, as t — > oo, 



F, '" l\ es dB s ) ^ U VH ^ 2H) N 



where N ~ Af(0, 1) is independent of B. 



Lemma 7 Let H > |. Then, as t — > oo, 



SB r e 6r — >0 mi 2 (ft), (14) 



t 

dse- 9s I dre 9r \s - r\ 2H - 2 — > 0. (15) 
Jo 

Proof of Lemma [BJ For any d > 1, si . . . Sd € [0, oo), we shall prove that, as t — > oo, 



o 



B Sl ,...,B Sd ,e- et e 9s dBA^ I B Sl , . . . ,B, d , v "T" y iV) (16) 



VHT(2H) 

e H 

which is enough to lead to the desired conclusion. Because the left-hand side in the previous 
convergence is a Gaussian vector (see proof of [Lemma 7, [5j]), to get (fl6|) it is sufficient to 
check the convergence of its covariance matrix. Let us first compute the limiting variance of 



9 



e~ et g e 0s dB s as t — > oo. We have 



E 



S 

Jo 



- et I e 0s dB, 



e 9s e er \ S ~r\ 2H - 2 drds 



= H(2H-l)e- 29t f f 
Jo Jo 

= 2H(2H - l)e- 29t f dse 9s [' dre 9r \s - v\ 2R - 2 
Jo Jo 

= 2H(2H - l)e- 29t f dse 29s f dre 
Jo Jo 

= 2H(2H-l)e- 29t j\re- 9r r 2H - 2 j\ 



-9r r 2H-2 



dse 



20.- 



H(2H 1) f f J2H-2 e -Br ^ T _ e -2Bt f J2H-2 Jr 



^ r 2H - 2 e 6r dr^J 



HT(2H) 

e 2H 



as t — > oo, 



because e~ 29t f* r 2H ' 2 e 0r dr < e~ et /* r 2H - 2 dr = (2 gf -> as * -> oo. 
Thus, 



lim £ 

£— i-oo 



HY(2H) 

e 2H ' 



Hence, to finish the proof it remains to check that, for all fixed s > 0, 

E{ B s x e- 9t ( 

t— >oo 



Urn E (b s x e- 9t e 9v dB,)j = 0. 



Indeed, for s < t, 



El B,x e- 0t / e 6v dB, 



= H(2H - l)e- 0t f dve 6v f du\ 
Jo Jo 

= H(2H - l)e- 0t f dve 9v f 
Jo Jo 

= H(2H - l)e- et f dve 6v f du\ 
Jo Jo 



2H-2 



du\u — v 



2H-2 



2H-2 



+ H(2H -l)e~ 0t dve 6v du{v - u) 
+ He- et f e^iv 2 "- 1 -(v- sf^dv 

J s 



2H-2 



It + Ji- 



lt's clear that I t — > as t — > oo. 

Using integration by parts, the therm J t can be written as 



Jt 








= He- 0t 






' S 

/ P ot 


= He- 0t 




< 





\2H-1\ 



\2if-li 



„2H-1 



2H- 1 



fe ev [{v-s) 

J s 



2H-2 „.2H-2 



v 2H ~ 2 ]dv 



2H-U + H(2H - l) e _et 



f ' e 0v (v-s) 2H - 2 dv 

J s 



A + Jl 
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Since H < 1, J} = fit 2 "- 1 - (t - s) 2 "- 1 } -> as t -> oo. 
On the other hand, 

J2 = g(2g~l) e - 9t f e 0v {v _ s) 2H-2 dv 



e 




H(2H- 


l)e es 


9 

H(2H- 


l)e es 


e 

H(2H- 


l)e es 


e 

H(2H- 


l)e es 



ri-s 

,-9t / „6u„.2H-2 



e / e u 



; „■•- , ~ 9t f e 9u u 2H - 2 du 

Jo 

e -0v {t _ v) 2H-2 dv 



du 



Fix it e (0, 1). The function i e [0, oo) i— >■ t 2H ~ 1 e~ 9tu attains its maximum at t = 2H e ~ l ■ Then 

sup^-^-^) = ce-^u 1 ' 211 < cu 1 - 211 , 
t>o 

with c = ( 2i ^~ 1 ) 2H ~ 1 . In addition. J** u 1_2ff (l - u) 2H ~ 2 du < oo, and for any u G (0, 1), 

f 2H-i e -et« (1 _ u) 2H-2 ^ o as t ^ OO. 

Therefore, using the dominated convergence theorem, we obtain that J 2 converges to as t — > oo. 
Thus, we deduce the desired conclusion. ■ 

Proof of Lemma [?1 Let us prove the convergence (|14|l . We have 



e~ 9t E / 8B s e' 9s / 5B r e 6r 



t r s 



r '"El I - ' SBrSB. 

2 



2 



2 



H 2 (2H-l) 2 e _ gt r e -e\v- s \ e -e\u-r\\ v _ u fH-^ s _ r fH-2 dudvdrds 
2 J\o,t] 4 



i[o,t]- 

= l - [E {B 2 )] 2 e-^ = \t^e-^ 
— > as t — > oo. 

For the convergence (|15j). we have 



H(2H-l)e-% f dse~ 9s f dre 9r 
Jo Jo 

< H(2H-l)e~% [ ds [ dr\s - r 
Jo Jo 



s - r\ 2H - 2 



2H-2 



+2H 

— > as t — ^ oo. 
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This finishes the proof. ■ 

Proof of the theorem [5J By combining © and Lemma @J we can write, 



e 6i 



(e t -e) 



e 0t f*dB s e s °f°dB r e- e r 
fields 

e^oo e-etf Q ee*dB s 

X 



e -26t J* e 26sQ ds ^ 

e~ et Jl5B s e- Ss J°5B r e Sr 
e- 26t /J e 26s £ 2 s ds 

-J?(2J? - 1 ^ ^ — 

e~ 2et /* e 2es C 2 rfs 

:= A e t xB e t -Cf-D e t . 
Using Lemme[5]and Lemma |31 we obtain that 

— >• 29 almost surely as t — > oo. 
According to Lemma [HI we deduce 



^ JHT(2H) N 



Moreover, 



law \ 

because , 9 ^ A/"(0, 1) and 2V ~ A/"(0, 1) are independent. 

Thus, by Slutsky's theorem, we conclude that 

x B| ^> 20C(1) as t — > oo. 
On the other hand, it follows from Lemma [3] and Lemma [7l that 

C? P -^>- as t — > oo, 

and 

£)f — !> almost surely as i — > oo. 
Finally, by combining the previous convergences, the proof of Theorem [5] is done. 
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